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It has been shown through the AdS/CFT correspondence that the energy loss of a 
fast quark in a strongly coupled = 4 SUSY Yang-Mills matter in the large N limit 
is given by the classical Lienard formula. I demonstrate that under quite natural 



Q^' assumptions about the dynamics of heavy ion collisions this leads to a universal (i.e. 

j-^' independent of the initial parton energy, but dependent on flavor and centrality) 

, upper bound on the energy of the partons escaping from the plasma. This bound is 

> 

00 ' a Yang-Mills analog of the Pomeranchuk bound in classical electrodynamics, where 

in ■ 

en ■ it is a consequence of radiation in a strong external field acting on a relativistic 
O 

\^ . charge. Since as a result the massive constituent partons are slowed down to a 

o , 

OO , velocity v < c, the angular distribution of the emitted radiation exhibits a broad 



"dead cone". If the properties of conformal and QCD matter at strong coupling 



^ ■ are qualitatively similar, the existence of this universal upper bound would have 

H ; 

y3 , dramatic implications for heavy ion experiments. 



One of the most striking results coming from RH. 
suppression of high transverse momentum particles 
from the parton energy loss has been predicted 5|, 
dense quark-gluon matter. The radiative energy loss 
magnitude of the observed suppression, see for examp 



3E 



heavy ion program is the observed 
3|, Such a suppression resulting 



7 as a signature of the formation of 
, Q] has been found to describe the 
le reviews . A definitive test of 
this mechanism can be performed with heavy quarks 1^ - since heavy quarks move with 



velocity v < c, the induced radiation must be depleted due to the "dead cone" effect - the 
vanishing of radiation intensity in the forward direction. This reduces the amount of energy 
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loss, and results in weaker suppression for charm and bottom quarks. The heavy-to-light 
ratios at high transverse momentum have thus been predicted to exceed unity [l^; this 
conclusion survives after the differences in the production mechanisms and fragmentation 
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functions for heavy and light quarks are taken into account 

It thus came as a surprise when RHIC experiments 13|, llJ] observed a strong suppression 
of the high transverse momentum electrons originating from the decays of charmed and 
beautiful hadrons. The magnitude of the observed suppression indicates that heavy and 
light quarks are attenuated very similarly in hot QCD matter, in sharp disagreement with 
the theoretical expectations. This remains true even in the range of transverse momenta 
that may be dominated by the decays of b quarks 13|, Moreover, the momenta of the 
heavy quarks also seem to be strongly deflected by the medium, as is evidenced by the 
observed elliptical flow (azimuthal anisotropy of the produced heavy quarks with respect to 
the reaction plane) 13, [isl, Q. 

While it is not yet entirely clear that a perturbative approach cannot be reconciled with 
the data, the heavy quark energy loss puzzle as well as some other observations in high 
transverse momenta phenomena IJ] provide an ample motivation to think about 

possible non-perturbative effects. Unfortunately, the set of tools that can be used to address 
the real-time dynamics of QCD in the strong coupling domain is limited. Therefore any 
information about the behavior of gauge theories at strong coupling is very valuable and 
may provide a hint on how to deal with non-perturbative effects in QCD matter. 

There has been a major breakthrough in the understanding of strong coupling dynamics 
of A/" = 4 SUSY Yang-Mills theory made possible by the AdS/CFT correspondence jl?!. [isl. 

. In particular, the strong coupling dynamics on the gauge theory side appears dual to 
the classical supergravity in AdS^ x 5*5 space. Of course, QCD and A/" = 4 SUSY Yang-Mills 
are very different theories, and conformal invariance of the latter is a crucial property which 
determines the metric of the AdS^ space. Conformal invariance results in the absence of 
confinement and asymptotic freedom in A/" = 4 SUSY Yang-Mills theory. Thermodynamical 
and transport properties of A/" = 4 SUSY Yang-Mills are also different; in particular the bulk 



viscosity (related to the scale anomaly of QCD |2l|, l22|, l23|) vanishes unless the conformal 
symmetry is broken 2d, l25|, U^■ Nevertheless one hopes that at least some aspects of 



the plasma behavior at strong coupling may be universal for A/" = 4 SUSY Yang-Mills 
and QCD at temperatures higher than the deconfinement and chiral symmetry restoration 
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temperature (but not much higher so that the couphng is still strong). Various properties 
of the plasma have been considered through the AdS/CFT correspondence, including the 



computation of the shear viscosity-to-entropy ratio in the strong coupling limit |27l . |28 |. 
The Langevin dynamics of drag force acting on a massive quark traversing the Af = 4 SUSY 
Yang-Mills plasma has been determined by solving the dual problem of string trailing in the 
AdSs Schwarzschild background |29|, [SO, l3l|; a comparison of AdS/CFT drag and pQCD 



predictions for observables in heavy ion collisions has been recently performed in [32|. A 
different approach proposed in |33;|] aims at matching Af = 4 SUSY Yang-Mills calculations 



to the perturbative QCD. 



including [34 



35, 
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37, 



'he related problem has been discussed in a number of papers, 



In an illuminating paper 



39j, it has been shown by Mikhailov that the energy lost by an 



accelerating external source in A/" = 4 SUSY Yang-Mills theory in the large limit is given 
by the following formula: 



Erad{T) 



2-K 



dt 



[v X a)' 



(1) 



where A = g'^N is 't Hooft coupling, v is velocity and a = -y is acceleration of the charge; we 
put c = 1. Remarkably, with the substitution 



27r 



(2) 



this is exactly the Lienard formula 



401 1 dating back to 1898 which describes the radiative 
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421. In the 



energy loss in external electromagnetic fields in classical electrodynamics 
case of electrodynamics, the linear dependence of the radiated energy on the path length is a 
consequence of the linearity of Maxwell equations: the emission of electromagnetic radiation 
at any time is not affected by the previously emitted radiation. 

Yang-Mills equations however are non-linear, and one expects, and indeed finds at weak 
coupling [3], a non-linear dependence of the radiated energy on the path length. This non- 
linear dependence is a purely non-Abelian effect: the amplitude of radiation at any given 
time is affected by the previously emitted radiation. Such non-linear, non-local dependence 



on the path traversed by the charge is also responsible 



br the evolution of jet structure 
well established and tested 



in the vacuum as described by DGLAP equations |43|, y 
experimentally. The linear local dependence of the energy lost by the quark in the strong 
coupling regime as given by ([1]) is thus a highly non-trivial result. 
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The method of 39| is based on considering a Wilson loop with boundaries given by the 
external quark and anti-quark sources; in AdS^ x 5*5 space, the quark and anti-quark are 
connected by a classical string with two boundaries. The worldsheet of this string is an 
extremal surface in AdS2 C AdS^, and the energy of the accelerating quark is lost to the 
excitations on this surface - non-linear waves. The propagation of the nonlinear wave on 
the string worldsheet in the large N limit is described by the classical sigma model 39|] , and 
the linear dependence of the final result ([1]) on the path stems from the integrability of this 



model 
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47|^. In the large limit the interactions with closed strings are suppressed, 



and so the entire lost energy can be attributed to the non-linear wave on the extremal 
surface stretched between the quark and anti-quark. It is yet unclear (at least to the present 
author) what excitation corresponds to this wave in the dual Minkowski space gauge theory 
language; nevertheless below we will attempt to give a qualitative picture based on the 
analogy with electrodynamics of strong fields. 

It should be noted that once the linear, local dependence on the path is established, 
Lorentz invariance and dimensional counting completely determine the structure of the for- 



mula for the lost energy ; 
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421] - it has to be proportional to the Lienard formula in classical 



electrodynamics. Indeed, the relativistic expression for the energy-momentum vector of 



the emitted radiation reads 



41 



I42 



with the substitution ([2]) 



2tt 



d^Xp d'^xP ^ 



(3) 



where r is the proper time; this is exactly Mikhailov's result 39|. This means that the 
coincidence of the result ([1]) with the Lienard formula cannot be considered as an evidence 
that the mechanism of energy loss at strong coupling is classical radiation; having this 
in mind, we will nevertheless for simplicity use the familiar language of electrodynamics 
and refer to the flow of the emitted energy and momentum described by ([T]) and ^ as 
"radiation". 

I — I [ — I 

The result of Mikhailov ^39|] has been verified and extended by Sin and Zahed 4^ and by 

Chernicoff and Guijosa Sin and Zahed in particular have argued that high momentum 
partons in a strongly coupled plasma would not be able to penetrate beyond the distance of 



^ If it appears that this hnear Abehan-hke formula holds in QCD at strong coupling, the reason may be 
the conjectured quasi- Abelian dominance proposed by 't Hooft |5c3|; this would imply a dynamical role 



for magnetic monopoles, also in the plasma 5ll. |52| 
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l/vrT 48|. Chernicoff and Guijosa have derived the expression for the dispersion relation of 



moving quark, and have considered also the case of finite temperature 



49|. 



In this paper I will show that under quite natural assumptions about the evolution of 
gauge fields in heavy ion collisions, the result ([T]) implies the existence of a universal (i.e. 
independent of the initial energy, but dependent on the mass of the parton and on centrality 
of the collision) upper bound on the final energy of the parton escaping from the strongly 
coupled matter. 

Let us begin by considering the special case of acceleration parallel to the velocity, a || tJ. 



Introducing 7 = \j\J\ — and the momentum p = 7771-1?, we get from ([I]) or (IHl) the 
expression for the power of radiation (radiated energy per unit time): 

dt ~ 271 -i^m^ Xdr) ' ^ ' 

Since dr = dt/'j, and dp/dr = 7 dp/dt = 7-F, where F is the force acting on the charge we 
get 

F'. (5) 



dt 2n m? 

The radiation power dErad/ dt in the case of a || t; is thus independent of the energy of the 
charge and is determined only by the magnitude of the external force, as is well-known from 
classical electrodynamics [4I, The situation when a \\ v is encountered for example 
when a quark jet propagates in the vacuum and is slowed down by the force of the string 
F = dp/dt = dE/dx = a, where a is the string tension. In this case the rate of quark energy 
loss as given by ([5]) does not depend on energy. 

The expressions found in 49] for the energy and momentum of the propagating quark at 
finite mass diverge as the value of the external force approaches 

2tt 



crit - 777^ (6) 



It is natural to interpret this in analogy with electrodynamics as the critical value of the 
force capable to produce quark-antiquark pairs from the vacuum. The limiting value of the 
field strength can be incorporated in a non-linear generalization of electrodynamics uniquely 
determined by Lorentz invariance and causality proposed by Born and Infeld 53|]. Indeed, 
on the string theory side the limiting value Fcru enters the Born-Infeld lagrangian on the D7 
brane; once F > Fcru, the creation of open strings becomes energetically favorable, and the 
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system becomes unstable [49j. The mass m in ([6]) should be thought of as the constituent 
quark mass related to the D7-brane parameter Zm = \/X/2TTm; the size of gluonic cloud 
around this constituent quark is Zm (see 49[ and references therein). 

As F — i> Fcrit, the energy loss of the quark jet according to ([5]) and (E]) is given by 



dt 



crit 



(7) 



and is due to the string fragmentation, i.e. creation of quark-antiquark pairs from the 
vacuum. The force acting on the quark thus arises from the polarization of the vacuum 
by the "supercritical" charge which is screened by the creation of quark-anti-qua rk p airs 
from the vacuum - this is the mechanism of quark confinement proposed by Gribov [54] ; for 



a review see 



55| . We thus have arrived at the interpretation of confinement force at zero 



temperature as being due to the energy loss of supercritically charged quarks in the vacuum; 
this interpretation provides a simple relation between the mass of the produced constituent 
quark m, the coupling A = g'^Nf. = X^-nag in Ai'c = 3 QCD, and the string tension a: 

„ 27r w? 



crit 



a. 



(8) 



Using m = 300 MeV and = 0.3 (corresponding to ^/X ~ 3.4), we get for the string tension 
cr ~ 0.85 GeV/fm - quite reasonable value consistent both with phenomenology and the 
lattice data. The corresponding size of the constituent quark is Zm = VX/2'Km ^ 0.35 fm. Of 
course. A/" = 4 SUSY Yang-Mills theory is not confining; but once confinement is introduced 
as an external force, the formula ([8]) can tell us whether this force is "critical" and would 
result in the production of constituent quark pairs. The numerical estimates performed 
above suggest that the confinement force is indeed "critical". 

The case of quark acceleration parallel to the velocity a \\ v applies to the fragmentation 
of the jets in vacuum, and to the energy loss in the direction parallel to the colliding beams 
in hadron collisions. However it does not apply to jets produced around mid-rapidity in 
nuclear collisions. High energy collisions are accompanied by the creation of strings, or 
longitudinal color fields, which would exert a force perpendicular to the velocity of the jet 
produced at mid-rapidity. Such longitudinal fields of "supercritical" strength (i.e. capable 
of producing^ quark-antiqua.r <: and gluon pairs) have been shown to emergejii heavy ion 



collisions 



56 



53, 



58 



59 



color glass condensate 



60|, l6l|, l78| from the saturated parton distributions 



62 , \6Wi in the 



64| . The produced longitudinal chromo-electric and chromo-magnetic 
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fields (termed "glasma" in j58j) possess non-zero Chern-Simons number [58l. l65l. l66l. 1671. |69| : 
the corresponding fluctuations of Cliern-Simons number fiave been measured in real-time 
Yang-Mills calculation on the lattice GSj. The presence of Chern-Simons number can induce 
in heavy ion collisions the violation of parity [70] (the possibility of spontaneous T and V 
violations has been considered in 71[; in the context of heavy ion collisions, it has also been 

been predicted to have a distinct signature [73i,y,Q 

- charge asymmetry with respect to the reaction plane, resulting in the electric dipole 
moment of the produced quark-gluon matter; see [TQ] for an overview. Recent prelimir iary 
experimental results indicate that this phenomenon may be present in RHIC data 77] . 
The propagation of charge in external classical flelds of the type considered above has been 
considered by Shuryak and Zahed 791], who have evaluated synchrotron-like radiation basing 
on the extension of electrodynamics treatment due to Schwinger [80] to strong coupling. 

The longitudinal color flelds will result in the acceleration perpendicular to the velocity, 
a A-V. In this case the formulae ([I]) or ([3]) yield for the radiation power 

dErad 1 



2 f?2 



(9) 



dt 271 m? 

This expression differs from ([5]) in one but very important way - it is proportional to the 



square of quark's energy E = 7m (note that this will be true for any flnite 



a and v). Eq(^ is of course well known in classical electrodynamics 4l|, |42(] where it 



a ng 



e between 



describes for example the energy loss due to synchrotron radiation. However in classical 



electrodynamics ([9]) has to be supplemented 
that the classical description still applies 



41 



jy the condition on the strength of the fleld so 
42( 1. Namely, the fleld strength Grest ("we use 



this notation to avoid the confusion with the force F) in the rest frame of the charge must 
obey the inequality 



m 



(10) 



In the frame where the charge moves with the velocity v, the fleld strength is G = Grest/li 
and the Lorentz force acting on the charge is F = eG. Therefore the condition ( JTOl) translates 
into the following condition on external force F: 



76^ 



or in terms of external fleld G 



7^ < 1- 
rri'^ 



(11) 



(12) 
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As emphasized in j41j, this condition does not prevent the ratio of the "radiation drag" force 
(Q to dominate over the external Lorentz force F; their ratio (we have substituted ([2]) for 
electrodynamics) is proportional to 

^I'F (13) 

and at large energies 'j ^ 1 will grow large even if ( fTTl) is satisfied. One is therefore justified 
to assume that the radiation drag force ([9]) is the dominant force acting on a relativistic 
particle at 7 ^ 1. 

Shuryak and Zahed [tqI have argued that in QCD the fields G are strong, the coupling 
e is large, and the external charge is massless, so the condition fll2p does not apply. They 
have thus concluded that Lienard formula ([T]) cannot be used in QCD [79]. However, the 
fact that equation ([T]) appears as the answer in the strong coupling relativistic problem 
in the AdS/CFT approach encourages us to take Lienard formula seriously. Indeed, the 
formula ([1]) holds in the ultra-relativistic lirnit; the coupling constant ^/X which replaces 



has been assumed large in the derivation 39(]; the mass m as discussed above has to be 



understood as a constituent mass; and the external force in AdS/CFT is limited only by 
the condition As emphasized above, Lienard formula relies only on locality, Lorentz 
invariance and dimensional analysis, and thus can be expected to have a wider range of 
validity than classical electrodynamics. 



Let us now come to the main point of this paper. As was noted by Pomeranchuk |41l. |81| . 
the formula (Q has a very interesting implication^. Indeed, the power of radiation is equal 
to the rate of energy loss, dErad/dt = —dE/dt = —dE/dx, where the last equality holds for 
a relativistic particle with f ~ c. We thus can re-write as 

dE y/X F^x 



dx 27r 



E', (14) 



where we have explicitly indicated the dependence of the external force F on the coordinate 
along the path. This differential equation can be easily solved by noting that dE/E"^ = 
—d{l/E), and integrating over the path of length L we get 



Ef Eo 271 Jo 



-^ = Tr + — / dx^^. (15) 



^ Pomeranchuk considered the radiative energy loss of cosmic ray electrons in magnetic field of Earth [81j . 
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As the initial energy of the quark Eq increases and goes to infinity, the final energy of the 
parton escaping from matter tends to a constant value! 

Replacing J dxF'^{x) by the product of the path length L and an average value F"^, we 
get for the upper bound on the energy of the parton 

27r 1 . . 

Ebound = l' ' 

Let us measure the magnitude of the external force in units z of the critical one given by 
(Q, F = zFcTiticai- We then get 

F (^^^ 
Abound 2n z^L' ^ ^ 

If we assume as before in our discussion of the string tension that the external force is 

critical, 2; 1, we get from f|T71) a formula which depends only the value of the coupling 

and the path length: 



Ebound — — r; (18) 

zn L 

where Ac is the critical value of the coupling corresponding to ([6]); this formula applies when 
the masses of the propagating quark and the quarks produced from the vacuum are the 
same, i.e. when the propagating parton is light. 

In a more general case we can measure the magnitude of the force acting on a quark of 
mass m in units Zf oi the critical value ^ for the creation of quarks of fiavor /, 

7t 

the formula (ITB]) then becomes 

f mV 1 1 

If) zj 



F = zf F/„j = z^ — m}; < z^ < 1, (19) 



Ebound — 7^ — -p- (20) 

2'K xnif zi L 



Let us now make some numerical estimates. Assume first that the magnitude of the 
external force F is given by the string tension, see ([8]). As before, we choose \/~\ ^ 3.4 
corresponding to ag — 0.3, and let 2; 1; we then get from (fTSll 



L(fm) ^ ^ 

this means that none of the light high transverse momentum partons would escape from the 
longitudinal string^. For charm quarks, taking m = 1.3 GeV, rrif = 0.3 GeV and — > 1 in 



^ Note that we are discussing the case of strong coupHng, and so perturbative processes are beyond the 
realm of this approach. 
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(EOj), we get 

35 

^ L(hnj ^^^^ 
Now let us take account of the fact that heavy ion colhsions can produce much stronger 
color fields than encoded in ([8]), as discussed above. It has been estimated that the magnitude 
of the produced fields at mid-rapidity is not much lower the critical one needed for the 



production of charm quarks 



82| . This leads us to take irif = rric, — 1 in fl20l) : we then 



get for charm the same estimate as we got above for light partons, f l2T|) 



^ ^^^^ 

and for beauty with ^4.5 GeV 

Hound - ^^f^) {Z^} 

We should of course admit that the estimates above are very rough, and depend crucially 
on the magnitude of the produced color fields. 

These estimates suggest that in strong color fields produced in relativistic heavy ion 
collisions at RHIC both light and charm quarks and gluons cannot escape from dense region 
of the produced matter. Of course this does not mean that there will be no high transverse 
momentum particles - they will be emitted from the surface of the produced fireball, leading 
to a universal normalized ratio of nuclear and proton-proton cross sections Raa for gluon, 
light and charm quarks almost independent of energy. A weak increase of Raa may result 
from the small amount of conventional absorption in the dilute "corona" surrounding the 
dense core of the plasma which is expected to go away at large transverse momentum in 
accord with factorization theorems of perturbative QCD. This leaves little room for observing 
the medium-induced modifications of the jet structure - either the jet is produced in the 
"corona" and is thus not modified at all, or it is produced in the dense core and is completely 
absorbed, with the final energy below the bounds given above. 

Based on the perturbative arguments, one is led to search for the jet modification in 
the central collisions of heavy nuclei. The bound presented here is inversely proportional 
to the length of traversed medium, and so suggests that the only hope to observe the jet 
modified by the medium-induced radiation is in peripheral collisions and/or in the collisions 
of lighter ions. The color fields, and thus the external force acting on the color charge, are 
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proportional to the saturation momentum squared, F Ql which grows with the centrahty 
as ~ ^pirt where Npart is the number of participant nucleons (see jsSi] for details). Since 
L ~ Np[l^ the bound may be expected to depend on centrality as Ebound ~ ^/ Npart- Of 
course, detailed calculations including the realistic nuclear geometry and taking account of 
the trigger bias effect enhancing the contribution of the smaller path lengths will have to be 
performed to get a rehable estimate of the bound. 

The color fields, and thus the external force acting on the color charge, are proportional 
to the saturation momentum squared F Ql which grows with energy the bounds will 
decrease at the LHC ener gy. According to the estimates of the saturation momentum at the 
LHC and RHIC (see e.g. [SJ]), the bound on the h quarks at the LHC will therefore change 
to 

2 

Abound - ^^f^) ^el/, {Z^} 

SO the suppression of h quarks at high should become a clearly visible effect. 

Let us note also that once the massive constituent partons are slowed down below the 
energy given by the bound, their velocity will become significantly smaller than the velocity 
of light. Therefore the emitted radiation must possess a "dead cone" of size 9 cone — m/E. 
The energy E of the radiating parton will decrease from the initial energy to the final energy 
given by the bound. To get an upper bound on the size of the dead cone we use E = Ebound 
and get 

m 2n ^ 
Ocone < -T. ^ ^ mL, (26) 

Abound V A 

where m is the constituent mass; we have used the expression (fT8l) for the bound on the light 
parton energy. Since the emitted energy-momentum is linear in the path of the parton, the 
resulting angular distribution will be a linear superposition of emissions from partons with 
energies ranging from some initial energy Eq to the final energy Ebound) with "dead cones" 
ranging from very small Ocone — ^/Eq to the large value fl26l) . m/EQ < 6 cone < ^/Ebound- 
Therefore the resulting dead cone can be quite broad; again, detailed calculations are needed 
to evaluate the shape of the angular distribution of the emitted radiation. 

To summarize, we have shown that if Lienard formula derived through the AdS/CFT 
correspondence holds in the strong coupling regime in QCD then under quite natural as- 
sumptions about the evolution of gauge fields in heavy ion collisions there should exist an 
upper bound on the final energy of the parton escaping from the strongly coupled matter. Of 
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course, the asymptotic freedom of QCD dictates that at some large transverse momentum 
the dynamics should become perturbative. However it is not yet clear at what scale the 
strong and weak coupling descriptions match in heavy ion collisions; an experimental study 
of the possible existence of the bound is needed to answer this question. 

This work was supported by the U.S. Department of Energy under Contract No. DE- 
AC02-98CH10886. 
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